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A new scheme is proposed which will permit electron spin resonance pulse 
techniques to be used to realize a quantum computer with a 100 qbits, or 
more. The computation is performed on effective pure states which correspond 
to off-diagonal blocks of the density matrix. Described is a scheme which 
very efficiently performs the preparation stage and which permits "pseudo- 
projective measurement" to be made on the output. With such measurements 
all members of the ensemble remain coherent. 

The demonstration by Shor (1) that a quantum computer can factor large numbers 
in a time which goes as a polynomial rather than the exponential of the number of digits 
has electrified interest in quantum computers. 

In a quantum calculation an initial state |0 >= |0 . . . > is acted upon by a unitary 
transformation U which reflects the calculation and which gives a result \R >= U\0 >. In 
order that the calculation make sense |0 > must be a perfectly specified state and it must 
be possible to determine \R>. The designer of a quantum computer is immediately faced 
by three non-trivial problems. First, since |0 > is some particular state, the entropy is 
zero and the system must be cooled to absolute zero. Second, the fundamental tenants of 
quantum mechanics insist a given quantum state \R > cannot be determined accurately 
by a single measurement. Repeated calculation cycles are necessary even before concerns 
of the signal to noise are addressed. Third, it is necessary to read the individual qbits. 

Cirac and Zoller (2) have suggested a physical realization of a quantum computer 
based on linear ion traps, however even supporters (3) of such devices as a proof-of-principle 
for quantum computing concede that these will never make a useful computer. Ingenious 
schemes for performing computations using standard liquid room temperature nuclear mag- 
netic resonance (NMR) techniques has been introduced by Gershenfeld and Chuang (4) 
and Cory et al (5). Again there is a serious scaling problem and it is difficult to imagine 
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a computer with much more than 10 qbits. As pointed out by Warren (6) in a ensemble 
of 10 23 computing molecules at room temperature there is essentially zero probability of 
finding a single example of any state designated to be |0 > for a 100 qbit register. 

An efficient ~ 100 qbit ensemble spin computer requires an almost complete polariza- 
tion of the spins and almost inevitably low temperatures, electronic spins, and the solid 
state. Commercial pulse electron spin resonance (ESR) spectrometers operating at ~ 100 
GHz are available. For spins with g = 2 this corresponds to a field of 3T. At a tempera- 
ture of IK, e hu o/k B T ^ iq-4 go f or a qbif. computer a large fraction of the computing 
molecules will in the ground state designated to be |0 >. A basic figure of merit Q = ujqt 
is the ratio of u>o the transition frequency (1/wo is the shortest switching time) and r the 
coherence time. As a general rule, light elements have the longest r, e.g., Li (7) or organic 
free (8) radicals have (7-factors very close to 2.00 and very long Ti ~ 10~ 3 s times, implying 
Q values as large as 10 9 in the solid state. For comparison DiVincenzo (9) gives Q = 10 7 
for NMR and Q = 10 13 for ion traps. However these latter figures are misleading since the 
implementation of quantum gates requires individual qbits to be addressed. Both linear 
traps and liquid NMR are very slow because Au;o the differences in frequencies are very 
small. In order to avoid this difficulty, it is envisaged here that each qbit comprises two 
spins. Rare earth, i.e, 4f ions such as Dy and Er often have g ~ 6 — 7 but somewhat shorter 
T\ times. It is imagined that, e.g., a free radical plus 4f ion bi-spin unit would be a basic 
building block with the 4f ion always in its very long lived ground state. The intra-unit 
exchange interaction between this 4f moment and the computing spin via differing bridge 
molecules would provide large shifts from g = 2. The inter- unit exchange provides the nec- 
essary interactions for qbit manipulation. The technology envisaged is that of molecular 
magnetism which has been much developed (10) over the past decade or so. 

The detailed realization of such a machine is not the principal concern here. The ques- 
tions asked assume that such physical systems can be engineered. Consider the following 
structure for a register of a quantum computer: There is a alternating chain of N comput- 
ing spins which have distinct resonance frequencies uj n even in the absence of interactions. 
To simplify the translation to the computational problem, for each computing spin the 
notation | |>= |1 > and | |>= |0 > will be adopted and it will be assumed that |0 > 
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corresponds to the ground state. Due to interactions, each spin sees its two neighbours, so 
that each resonance at u n is split into four addressable (well separated) lines with frequen- 
cies uj® , uj® 1 , iv^ and uj^ 1 as the neighbours are 00, 01, 10 and 11 respectively. Here 10 
means that the spin n — 1 is 1 while the n + 1 spin is 0. Because of the alternation of the 
interaction u® 1 ^ u^ . A it pulse at any one of these frequencies amounts to the execution 
of a generalized Toffoli gate, i.e., such a pulse, e.g., at the frequency ui^ will flip the nth 
spin if, and only if, its neighbors are "1" and "0". This gate along with arbitrary spin 
rotations can be effected by standard pulse resonance technique and permit all quantum 
gates to be constructed. This has been adequately dealt with by others (4,5,11,12) and 
will not be described here. 

Even with the considerable advantage of low temperature electronic spin systems there 
remain two very important problems in using finite temperature ensembles for quantum 
computing. First, even if they are relatively few, it is still necessary to eliminate the 
signals from molecules not in the ground state |0 >. The result \R >= U\0 > will be 
a linear superposition of many states and the weight of a given state in \R > will very 
often be much smaller than the weight of unwanted signals from other members of the 
thermodynamic ensemble. Necessary is a preparation (or final averaging) procedure in 
order that the observable signal originates only from the computation pure state |0 >. Most 
of the procedures suggested to date (4,5,12,13) are for small high temperature ensembles. 
They typically become exponentially long for large registers negating the advantage of 
performing a quantum computation. Potentially non-exponential schemes do exist (13) 
although no specific implementation has been suggested for a low temperature ensemble. 
The principle of these earlier techniques is quite different from that proposed here. The 
usual idea is to reshape the entire density matrix to be an effective pure state of the form 
(1/N) +p\0 >< 0|, while here the aim is to create off-diagonal blocks \a > |0 >< 0| < — cr|, 
a = 0, 1 which correspond to an ESR signature from a "signal" spin reflected by a. The 
second problem has to do with projective measurements. Quantum computations typically 
involve measuring the state of \R > and, after the principles of quantum mechanics, this 
involves a random projection. If such measurements were possible using NMR, or ESR, 
pulse techniques this would leave different members of the ensemble in randomly different 
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states. A remedy for this difficulty has been suggested (4), however it remains the case that 
such perfectly projective measurements do not exist in the repertoire of NMR, or ESR, 
pulse techniques, although a clever such scheme does exist for ion traps (14). It is always 
possible to imagine deterministic equivalent of a given quantum algorithm although no 
explicit such scheme has yet to be presented of the factoring problem. Here is proposed a 
general method by which to perform "pseudo-projective" measurements and which permit 
the existing algorithms to used directly. 

The basic idea of the new preparation scheme is to perform a simple quantum com- 
putation on all the states of the ensemble which permits the free induction ESR resonance 
signal to be reduced to that coming from the state |0 >. A signal spin, the ESR of which 
is monitored, is added as the first spin of a AT + 1 spin register, thus if \A > is an arbitrary 
(N spin) state, the state of the whole register is \a > \A>, a = 0,1. The equilibrium en- 
semble gives no (free induction decay) ESR signal. A signal corresponding the transitions 
|<7 ><^> | — a >, for states of the form \a > \0A >, is turned "on" by ir/2 pulse at u)®. The 
standard "spin echo" technique comprises a sequence of gradient- 7r-gradient pulses (15). 
A field gradient applied to the sample causes the spins on different molecules to precess at 
different rates and this destroys the ESR signal, however the signal is recovered following a 
time conjugating n pulse if an identical gradient pulse is applied to re-focus the spins. The 
principle of the preparation stage is, between gradient pulses, to perform a computation 
on all the states of the ensemble which time conjugates the signal spin if, and only if, the 
state | A >= |0 >. Conceptually the simplest scheme would be to consider the signal spin 
as a separate single bit register s± and perform a quantum calculation equivalent to the 
pseudo code 

If A = 0, s 1 = s 1 + l (mod 2) 

The result of a computation is |cr > |i2>=|cr>C/|0> where U is the relevant 
unitary transformation which has no effect on the signal spin. Projective measurements 
are replaced by what might be called "pseudo-projective" reversible equivalents. The 
principle (again the practice is somewhat different) is similar to the preparation stage. 
The presence of a state \a > \S >= \a > |a 2 , cr 3 . . . > is detected by a signal spin ESR 
by performing the equivalent of a calculation on all states \a > \A >, between two gradient 
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pulses, which reverses the signal spin if, and only if, A = Sq, i.e., the pseudo code is: 
If A = S , s 1 = s 1 + l (mod 2) 

The same pulse sequences in the opposite order return the system to the state before 
the interrogation process and hence this reading process is reversible. Usually the result 
cannot be determined by a single such measurement and measuring every amplitude takes 
an exponentially long time. Described here are routines which take, at the most, ~ iV 2 
pulses to read the result for either the Shor or Grover (16) algorithms. 

While the advantage of a reversible "read" process are considerable, there is also a 
major potential disadvantage over projective measurements. After a perfect projective 
measurement (if such a thing really exists), of, e.g., spin orientations corresponding to the 
set of operators <r*, n = 1, N + 1, the system finds itself in an eigenstate n^i 1 \ a n > °f 
these operators. The weight of this state in original wavefunction is reflected not in the 
strength of the associated "signal" but rather by the probability of finding the particular 
signature corresponding to the set {u^}. In a straightforward reversible measurement of a 
quantum state the strength of the signal must be proportional to the weight squared of the 
particular state. In a N + 1 qbit register, and the worst case when all states have roughly 
the same weight, the strength of a reversible measurement is ~ 2~ N as is the probability 
of a particular result of a projective measurement. An efficient quantum algorithm is 
explicitly constructed to have a \R > with much less quantum entanglement, however this 
can remain a difficulty for reversible reads. In this regard Grover's algorithm presents no 
problem since the answer dominates the register to be measured. The Shor scheme requires 
more discussion which will be taken up again below. 

The first step is to describe an efficient scheme which reduces the ESR response of 
the signal spin to that coming from the two computational \a > |0 > states. This uses no 
ancilla spins other than the signal spin. As described above, the signal comes only from 
states of the form \a > \0A >. After a first gradient pulse, are 7r pulses at the frequencies 
a>*°; n = 2, N. A second less selective it pulses sequence, which only "looks left", has 
frequencies n = N, 2 where X = Oorl. (It is implied that pulses are applied at, or 
cover, both and cu^°.) Last is a n pulse at u\. For the computational states the net 
effect is |1 > |0...00 >-> |0 > |1...10 > and |0 > |0 . . .00 >-> |1 > |1 . . . 10 >. If the last 
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spin was reversed this would be a complete time conjugation, however this is not necessary 
since both states de-phase in the same manner for a given molecule and so a second 
identical gradient pulse does re-focus the spins and the ESR of the signal spin is recovered. 
For the non-computational states there are two cases. If the first "1" occurs on the third 
spin, the net effect is |1 > |0E4 >-> |1 > \00A" > and |0 > |0E4 >-> |1 > \10A" > 
where A" is some unique mapping of A. Generally the first "1" occurs after M zeros 
|1 > |0...00E4 >^ |0 > |1...KXL4" > while |0 > |0...001A >^ |1 > |1...110A" > 
where . . . replaces a series of digits "1" . In all cases the resulting states differ at the qbit 
positioned two places before the "A" and as a result the final states dephase differently 
during the second gradient pulse and the ESR signal is not recovered. In addition, for the 
first case, the signal spin is not properly conjugated. 

This preparation stage permits the extraction of the result of the quantum compu- 
tation on the pure states \a > U\0 >. The transformation U does not change the signal 
spin ESR. The weight of any given state |a"2 • • - o~n > in U\0 > might be determined by 
another gradient-7r-gradient sequence, however it would take an exponentially large time 
~ 2^ to measure the amplitude of every state. There are probably as many strategies for 
rendering practical such a scheme as there are quantum algorithms (so finding them is not 
a big task as of writing). 

The Shor routine searches for the period r of f(a) = x a ModiV where N is a number to 
be factored and x < N is randomly chosen. The result is \R >= (1/ 'y/r) J2p=o \pT > \F p >, 
and the sought for r = w/T where w corresponds to the size of the registers and where T 
is the period, i.e., the smallest non-trivial number in the first register. The (normalized) 
\F p > need not be read. This smallest number has the largest quantity, N, of leading zeros, 
i.e., it is of the form \a > \0 . . .OA >= \si = a > \Sp = . . . > Q \A >p. The pseudo-code 
which conjugates the signal spin for such a state is 

If S/3 = 0, s 1 = s 1 + l (mod 2) 

This program line is run with N = 1, 2 . . . until it fails, with iV = iVo, to conjugate 
the signal spin. The last such pulse sequence is then run backwards to recover the result 
before the failure. The smallest number then begins \a > |0. . .01 A' > with N — 1 zeros 
and the search is continued for the smallest A' : etc. 
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Again a more efficient equivalent process is possible. The \a > \0A > will exist if a 
significant signal at uj® is observed. Assuming this is the case, a finite amplitude for the 
states | a > |00A > is sought. First the signal from the \1A > states is "killed" by a n/2 
pulse at u\. The remaining task is relatively simple since it is known where the first "error" 
will occur. As above the process starts with a gradient pulse, followed by ix pulses at uj^ 1 
and u^. The net changes are |1 > |01A >^ |1 > \l\A >, |1 > |0(L4 >^ |0 > |0(L4 >, 
|0 > |0(L4 >^ |1 > |0(L4 > and |0 > |01A >^ |1 > |01A >, i.e., this performs the 
necessary conjugation on the signal spin. The final gradient pulse will restore the ESR 
signal from the states \a > \00A > while leaving the states originating from \a > \01A > 
with the signal turned off. 

The generalization to \a > \0...0A >, where there are A zeros, contains the n 
sequences n = 2, A — 1; uj} followed by cu^ , n = A — 1, 2 which resets the state 
|1 > |0 . . . 00A > while having the effect |1 > |0 . . . 01A >^ |1 > |11 . . . 11A > but leaving 
the states |0 > |0...0CL4 > and |0 > |00...01A > unchanged. Clearly the sandwiching 
gradient pulses destroy all but the signal from \a > |0 . . . OA >. The generalization to the 
search for other specific states is evident. 

The signal corresponding to {l/^/r)\T > is reduced by a factor ~ 1/r where r < A, 
A being the number to be factorised. It is possible to detect ~ 10 5 equivalent spins, so 
if r > 10 18 even for a 10 23 element ensemble the signal for a ~ 100 qbit register will 
become unreadable. However a pseudo-collapse of the wavefunction can be accomplished 
by changing, a little, Shor's algorithm. The usual calculation, U, causes \a > \b >^ \a > 
\x a b >, where again a; is a randomly chosen c-number less than A. The pseudo-collapse is 
accomplished by modifying the second register so that the effect olU is \a > [|1 > +|x > 
+ \x 2 >+... + \x n >+... + \x s >] -> \a > [\x a > +\x a+1 > +\x a+2 > + . . . + \x a+n > 
+ . . . + \x a+s >] where s < r. This modification can be performed efficiently but requires 
an extra register. Following an initial failure to read \T > with s = 0, a second calculation 
with s ~ 10 18 must succeed. 

For the Grover algorithm the answer dominates U\0 >. First the signal spin is ob- 
served. There will be a single strong line which determines the leading digit in U\0 >. 
Imagine, e.g., that this implies U\0 >~ \1A >, using the method of the previous few para- 
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graphs the amplitude of \11A > is determined. If, e.g., this is small then the state begins 
\10A > and the last interrogation sequence is run backwards to recover the amplitude for 
this state and then |101A > sought for, etc. (An alternative but somewhat more "rough 
and ready" method is to apply a small angle, say 7r/20, pulse to all spins. This will turn 
"on" an ESR signal which corresponds to the desired state, e.g., if U\0 >~ 1 10100010 . . . > 
then the n = 2 spin will give a free induction signal at the two positions uj* 1 and then at 

, ,11 , ,00 , ,10 \ 

Error correction represents an interesting challenge. In order to remove the entropy 
generated by errors, such schemes (17) use projective measurements in an essential fashion. 
Replacing these by the pseudo-projective equivalents would cause the signal strength to fall 
off exponentially. However it is not hard to find alternative means by which to perform the 
correction using conditional unitary transformations. The entropy must be carried away 
by ancillary qbits and then removed by cooling. A full discussion of these possibilities is 
too lengthy to be reproduced here. 

The author thanks E. Knill for helpful comments. 
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